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A sufficient condition for Gaussian departure in turbulence 
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The interaction of two isotropic turbulent fields of equal integral scale but different kinetic energy 
generates the simplest kind of inhomogeneous turbulent field. In this paper we present a numerical 
experiment where two time decaying isotropic fields of kinetic energies Jf?i and E2 initially match over 
a narrow region. Within this region the kinetic energy varies as a hyperbolic tangent. The following 
temporal evolution produces a shearless mixing. The anisotropy and intermittency of velocity and 
velocity derivative statistics is observed. In particular the asymptotic behavior in time and as a 
function of the energy ratio E1/E2 — > 00 is discussed. This limit corresponds to the maximum 
observable turbulent energy gradient for a given E\ and is obtained through the limit £2^0. A 
field with E1/E2 — > 00 represents a mixing which could be observed near a surface subject to a 
very small velocity gradient separating two turbulent fields, one of which is nearly quiescent. In 
this condition the turbulent penetration is maximum and reaches a value equal to 1.2 times the 
nominal mixing layer width. The experiment shows that the presence of a turbulent energy gradient 
is sufficient for the appearance of intermittency and that during the mixing process the pressure 
transport is not negligible with respect to the turbulent velocity transport. These findings may 
open the way to the hypothesis that the presence of a gradient of turbulent energy is the minimal 
requirement for Gaussian departure in turbulence. 

PACS numbers: 47.27+, 47.51+ 



I. INTRODUCTION 

A turbulent shearless mixing layer is generated by 
the interaction of two homogeneous isotropic turbulent 
(HIT) fields, sec definition diagrams in figures [T] and 
[2] and the flow visualizations in figure [3j This kind of 
mixing is characterized by the absence of a mean shear, 
so that there is no production of turbulent kinetic en- 
ergy and no mean convective transport. The turbulence 
spreading is caused only by the fluctuating pressure and 
velocity fields. The inhomogeneous statistics are typi- 
cally due to the presence of the gradients of turbulent ki- 
netic energy and integral scale. The shearless turbulence 
mixingwas first experimentally investigated by Gilbert 
(1980)13 and by Veeravalli and Warhaft (1989) by 
means of passive grid generated turbulence. Later on, 
numerical investigations were carried out by Briggs et &1. 
(1996) and Knaepen et ai. (2004) 0, and more re- 
cently by Tordella and Iovieno (2006, 2007) 0,0. All 
these studies considered a decaying turbulent mixing. 

In all studies, apart from that of Gilbert, where the 
turbulent energy ratio was very low, the mixing layer 
was observed to be highly intermittent and the trans- 
verse velocity fluctuations seen to have large skewness. 
Across the mixing the distributions of the second, third 
and fourth order moments collapse when the mixing layer 
width is used as lengthscale 0,111, • 

In passive grid laboratory experiments the gradients of 
integral scale and kinetic energy are intrinsically linked. 
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In past studies the ratio of the integral scale of the inter- 
acting turbulence fields was in the range 1.3 - 4.3 
with a ratio of kinetic energies in the range 1.5 - 23 
. In numerical or active grid experiments these two 
parameters can be independently varied. 

In the present study, a mixing configuration in which 
the integral scale is homogeneous is considered. The ra- 
tio of the turbulent kinetic energies has been chosen as 
the sole control parameter and is varied from 1.5 to 10 6 , 
Re\ of the high turbulent energy field is 45. The aim of 
this study is to show the intermittent behavior of such 
a configuration that in the past was considered to have 
almost Gaussian velocity statistics. This interpretation 
was motivated by the absence of both a kinetic energy 
production and an integral scale variation, two typical 
sources of intermittency and was also supported by lab- 
oratory observations carried out in the absence of a suf- 
ficiently high kinetic energy gradient 0. Another aim 
of this numerical experiment is to reach the asymptotic 
condition where the kinetic energy ratio £ = E\ / E2 goes 
to infinity. This last condition is relevant in applications 
concerning the diffusion of a turbulent field in a region 
of quiescent fluid, where extreme bursts of rate of strain 
and vorticity can be expected 0. The presence of such 
events is shown by high values of skewness and kurtosis. 

A description of the numerical experiment is given in 
section II. Data on the degree of anisotropy observed 
in the second and third order velocity moments are de- 
scribed in section III, where an interpretation based on 
Yoshizawa's hypothesis is also given. In section IV we 
present the two types of asymptotics considered: the tem- 
poral asymptotics of the second and third order velocity 
moments, and the asymptotics with respect to the turbu- 
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responding to the high energy field, is equal to 45 for both 
the spatial discretization of the direct numerical simula- 
tions (DNS). 

In the initial condition, the two isotropic turbulent 
fields are matched by means of a hyperbolic tangent func- 
tion. This transition layer represents 1 /40 of the 4ir do- 
main, and 1/80 of the 87r domain. The matched field 
is 



u(x) = ui(x)p(x) + u 2 (x)(l -p(x)) 



(1) 



FIG. 1: Scheme of the flow. Direction x is the mixing di- 
rection. The high energy (Ei) and low energy (E2) regions 
are separated by mixing layers of conventional thickness A(t) 
defined by mapping the low energy side of the mixing layer 
to zero and the high energy side to one. A(t) is equal to the 
distance between the points with normalized energy values 
0.25 and 0.75 0,0. 
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FIG. 2: Scheme of the flow. Reference frame: 2/1,2/2 are nor- 
mal to x, that is the direction of the flow inhomogeneity. The 
flow is homogenous in all planes normal to this direction. 



lent kinetic energy ratio of the velocity skewness, mixing 
penetration and kurtosis. In addition, a smaller set of 
data on the temporal asymptotics of third and fourth or- 
der moments of the velocity derivative is also discussed in 
this last section. The concluding remarks are presented 
in section V. 



II. NUMERICAL EXPERIMENT 

Navier-Stokes equations are numerically solved with a 
fully dealiased (3/2 - rule) Fouricr-Galcrkin pseudospec- 
tral method The computational domain is a par- 
allelepiped with periodic boundary conditions in all di- 
rections, see fig.l. Tests were performed on a 47r(27r) 2 
parallelepiped domain with 256 x 128 2 points. Further 
tests with a 8ir(2ir) 2 parallelepiped with 512 x 128 2 points 
were used to obtain an estimate of the numerical accu- 
racy. The Taylor-microscale Reynolds number Re\, cor- 
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where the suffixes 1, 2 indicate high and low energy sides 
of the mixing respectively, x is the inhomogencous direc- 
tion, L is the width of the computational domain in the 
x direction. Constant a in (??) determines the initial 
mixing layer thickness A, conventionally defined as the 
distance between the points with normalized energy val- 
ues 0.25 and 0.75 when the low energy side is mapped to 
zero and the high energy side to one. When a = 12tt the 
ratio A/L is about 0.026, for a = 20ir the ratio A/L is 
about 0.015. These values have been chosen so that this 
initial thickness is large enough to be resolved but small 
enough to have large regions of homogeneous turbulence 
during the simulations. This technique of generating the 
transition layer is analogous to that used in Briggs et 
al. (1996)0, and Knaepen et al. (2004)0). The match- 
ing on which the initial condition is built up is a linear 
superposition of the two isotropic fields as indicated in 
equations (jTJ) and (??). A set of statical preperties of the 
high kinetic energy HIT field is shown in Table 1. Since 
the low energy field U2 is obtained by multiplying the 
initial velocity field u± by a constant, the numerical ex- 
periment carried out by mixing these fields is a turbulent 
mixing with different energies but of equal integral scale. 
It should be noted that, by doing so, the mean pressure 
along the mixing direction is not constant, However, the 
mean pressure gradient is opposite to the gradient of tur- 
bulent kinetic energy and thus no mean velocity field is 
generated, see the appendix. Examples of the shearless 
mixing obtained in this way for direct numerical simula- 
tion can be found in Q and Q. The initial spectra of 
the two HIT fields are shown in figure 2J In this figure 
the temporal decay of the two isotropic turbulent fields 
is shown together with, as a reference, the decay of the 
homogeous and isotropic turbulence simulated in one of 
the computational domain used to simulate the turbulent 
shearless mixing ((27r) 2 x 87t, 128 2 x 512). In figure[3]the 
estimate of the time instant where the self-similar decay 
of the mixing starts is also shown. 

Let us now consider the flow symmetry. It can be seen 
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FIG. 3: (Color online) Visualization at two time instants of contours of kinetic energy E(x,yi, 7/2, t)/Ei(0) in a plane at constant 
2/2, Ei/Ez = 6.7, Re x = 45: (a) t/r = 0.8, (b) t/r = 2.5. 

that a shcaiicss mixing is a flow in which only one direc- 
tion of inhomogeneity is present, as a consequence any 
plane normal to the inhomogcncous direction is homoge- 
neous. This corresponds to a cylindrical symmetry. See 
the reference frame scheme in figure O 

The time integration is carried out by means of a four- 
stage fourth-order explicit Runge-Kutta scheme. Statis- 
tics are obtained by averaging over planes normal to the 
inhomogeneous direction, sec figure 

The initial conditions were generated from the homo- 
geneous and isotropic turbulent field produced by Wray 
in 1998 which is a classic data set often used in liter- 
ature. 

A posteriori, it is possible to obtain numerical accuracy 
estimates. The raw data by Wray has an inhomogeneity 
level on the kinetic energy of about ±8% and skewness 



and kurtosis values slightly different from those of the 
statistical equilibrium (0.02 ±0.12 instead of and 2.8 ± 
0.2 instead of 3, respectively). As far as our set of direct 
numerical simulations is concerned, the increase in width 
of the computational domain from 47r to 87r (from 256 
to 512 grid points) allowed an estimate of the relative 
accuracy to be obtained. For the maximum values of the 
distributions across the mixing, the accuracy is of about 
5% for the skewness, and of about 8% for the the kurtosis. 

In figure [51 which summarizes the results regarding 
the maximum values reached by the velocity skewness 
and kurtosis within the mixing and the results about the 
penetration, it can be seen that the simulations with ini- 
tial A/L = 1/40 and 1/80 yield data which collapse in 
a satisfactory way. On checking the symmetry of the 
numerical solutions, which, due to the periodicity of the 
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FIG. 4: Turbulent kinetic energy decay of the two interacting 
isotropic turbulent fields (Ei high energy, E2 low energy) at 
re\ = 45 and, in the inset, the corresponding initial energy 
spectrum. Data from a homogeneous and isotropic turbulence 
(E) simulated in a (27r) 2 x 8n domain (128 2 x 512) have been 
shown for comparison. The initial spectrum is equal to the 
spectrum of the high energy region in the mixing. 
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FIG. 5: Anisotropy of the turbulent second and third order 
moments at the centre of the mixing layer. The horizontal 
dashed line in part (a) indicates the isotropic reference value, 
the horizontal dotted line in part (b) indicates the estimate 
of the asymptote value. 



Sdu/dx Kdu/dx Sdu/dy,l ^du/dyS 

-0.42 ±0.08 3.61 ±0.2 -0.40 ± 0.08 3.53 ± 0.2 



TABLE I: Statistical properties of the high energy HIT field 
ui at t — 0. Ei is the normalized turbulent kinetic energy, 
Si and K\ are the velocity skewness and kurtosis, Sg u /g x and 
Kdu/dx are the velocity longitudinal derivative skewness and 
kurtosis, Sg u /g y ^ and Kg u /g Vt i are the velocity transversal 
derivative skewness and kurtosis. The field iti was obtained 
from the data base by Wray 1998 0. Note, that these statis- 
tical properties are the same for all the considered low energy 
fields U2, since they were obtained by multiplying the initial 
high energy field by a constant. 

boundary conditions, contain two mixings, see scheme in 
figure[l] it was verified that the doubling of the computa- 
tional domain induces a decrease of the asymmetry from 
10% to 5% for the skewness and from 20% to 15% for the 
kurtosis. 



III. ANISOTROPY AND YOSHIZAWA'S 
HYPOTHESIS 

In isotropic turbulence the normalized second order 
moment of the velocity components, normalized with the 
sum u 2 + v\ + t>2, is 1/3, whilst the third order moment 



is zero. In the present flow the field anisotropy develops 
during the mixing process. The value of the normalized 
moments vary in time and reach an asymptotic value 
after few time units, sec figure [5] The time unit r is 

1 /2 

defined as t = £(0) / (0)), where I is the integral scale, 
here uniform across the mixing, and E\ is the turbulent 
kinetic energy of the high energy side of the mixing. 

An initial turbulent energy gradient VE = {E\ — 
£"2)/(2A) corresponds to each value of £ = E\jE2. The 
width A is defined by mapping the low energy side of the 
mixing layer to zero and the high energy side to one, and 
it is equal to the distance between the points with en- 
ergy values 0.25 and 0.75, as in the paper by Veeravalli 
& Warhaft (in the following referred to as V&W). 
The turbulent energy gradients can be normalized by the 
value of the high energy field, and by the value of the mix- 
ing thickness A(i). It should be noticed that by doing 
so, the normalized gradient value has the upper limit of 
0.5, which is reached in the limit for E2 going to zero. 

In figure [HJa) the time evolution inside the mixing of 
the second order moment u 2 /(u 2 + v± + v 2 ) is shown. 
After a linear growth the curves bend toward the asymp- 
totic value, which is in the range 0.37-0.39 for a kinetic 
energy ratio growing from 4 to 10 4 (this corresponds to a 
normalized gradient of turbulent kinetic energy from 0.37 
to 0.50, or, by supposing a mixing in air with a Re\ = 45 
in the high energy side, to a dimensional gradient from 
1.8 to 2.4 m/s 2 ). 
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As a consequence of the cylindrical symmetry of 
this mixing, it follows that the second order moments 
v\j(u 2 + v 2 + v 2 ), v 2 j{u 2 + v\ + v 2 ) are equal and range 
from 0.315 to 0.305 when £ varies from 4 to 10 4 . The 
anisotropy level, defined as the difference between the 
second-moment values referred to the isotropic value, can 
be considered mild (16 % for £ = 12, 25 % for £ = 10 4 ) 
given that the accuracy in the original data base used 
to build the initial condition [TJ] is of about 8% as far 
as both the homogeneity and isotropy are concerned. It 
should be considered that this level of initial accuracy 
of homogeneity and isotropy is excellent in nominal HIT 
numerical fields. In higher resolution fields (1024 3 ) the 
accuracy is analogous [10| . 

Figure [Ha) indicates that the value 0.39 for u 2 /{u 2 + 
vf + v%) is reached by increasing £ from 12 to 10 4 . 
This value can be considered as an approximation of the 
asymptotic value attainable by increasing the turbulent 
energy gradient. 

It is important to note that in literature concerning the 
shearless mixing, almost all authors report a near homo- 
geneity in the second-order velocity moments regardless 
the observation method used, numerical or laboratory 

an. 

The anisotropy of the third-order velocity moments is 
more enhanced than that of the second-moments. This 
can be observed in figure EJb), where the time evolution 
of the third order moment u 3 normalized wit h the tota l 
kinetic energy flow in the mixing direction, u 3 +vfu+v 2 u, 
is plotted. The estimate of the temporal asymptotic value 
we obtained is 0.53 ± 0.03 and does not depend on £. If 
the level of anisotropy is defined as the difference between 
the third moments divided by their mean, an anisotropy 
of 80% is obtained. This means that, for all the energy 
ratios, nearly one half of the turbulent kinetic energy flow 
across the mixing is due to the self transport of u 2 . Let us 
note that at the initial instant, when the mixing process 
starts, the quantity u 3 /(u 3 + v\u + v\u) is not defined 
because both the numerator and the denominator are not 
defined. This is numerically verified through the large 
dispersion of the initial values associated to different £ . 
Of course, this dispersion is also due to the non perfect 
homogeneity of the HIT data base used to build the initial 
condition, see section II. The data dispersion is however 
reduced as the mixing process advances. After 6 times 
scales is less than 10%. 

It is possible to analyze this result by means of sim- 
plifying hypotheses currently found in literature - (a) 
the pressure transport is almost proportional to the con- 
vective transport associated to the fluctuations (Lumley 
1978[n|, Yoshizawa, 1982, 2002, El, [3), - (b) the dissi- 
pative scales are nearly isotropic [14| , and - (c) the second 
order moments are almost isotropic as observed in shear- 
less turbulent mixings and also confirmed by the present 
numerical experiment, as discussed above. 

Let us now consider the one point second order moment 
equations 



d t u 2 + d x u 3 = —2p d x pu + 2p pd x u — 2e u 



d t v 2 + d x v 2 u = 2p- 1 pd yi v l - 2e Vi + ud£v$, i = 1, 2 (4) 

where u is the fluctuating velocity in the inhomogeneous 
direction x, Vi,V2 are the fluctuation components in the 
plane normal to x and e u ,e Vi are the dissipation terms 
in the mixing and normal directions, respectively. 



The pressure strain terms pd x u and pd Vi Vi, in the ab- 
sence of a mean flow, are of the order of j (pu jUj — | pbSij ) , 
see for instance Monin & Yaglom, 1971 [15| (Volume 1, 
equation 6.12, page 379), where e is the total dissipation 
and b is the turbulent kinetic energy per unit of mass. 
Since, as previously explained, experiments show no ap- 
preciable difference in the second order moments in the 
mixing, see condition (c) above, the pressure strain terms 
are neglected. 

Condition (a) implies that we can write 



pu = ap- 



u 3 + 2v 2 u 



(5) 



for any value of position x along the mixing and for any 
time instant t. The difference between equation ([3]) and 
equation (0| gives 

d t {^-vj)+d x {^-vfu) « -2p-'d x pu-2(e u -e Vi ). (6) 

By condition (b) u 2 w vf and by condition (c) £ u ~ e Vi . 
Thus, the unsteady term on the left hand side as well as 
the second term on the right hand side can be neglected 
and it follows that 

d x (tf-vfu) ~-2p- x d xV u- (7) 
Integration of ([7]) with respect to x leads to 

u 3 — vfu ~ — 2p~ 1 pu + C, 

but, considering that all quantities in this equation van- 
ish outside the mixing (i.e. for x — > ±00), the integration 
constant C is equal to zero. Thus 



u 3 — vfu ~ — 2p 1 pu, 



(8) 



By inserting the previous relation into (O, it is possible 
to write 



vfu = f3u 3 , (3 



1-a 
1 + 2a 



(9) 



Then, by defining $ the proportion of the turbulent ki- 
netic energy flow associated to the u fluctuation, it fol- 
lows that 



$ = — 



1 



[u 3 + 2vju) 1 + 



(10) 
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We have computed the constant a for the present ex- 
periments and found, in asymptotic temporal condition 
and for £ e [12, 10 4 ], an average value of 0.37 ±0.03. This 
gives j3 ~ 0.36 and $ ~ 0.58 This last value contrasts 
with our numerical experimental value of $ = 0.53 ±0.03 
shown in fig. 5 b. 

We have verified that a and (3 remain almost constant 
during the decay and when varying the shearless mix- 
ing parameter £, a fact which confirms that the pres- 
sure transport correlation is almost proportional to the 
convective transport associated to the fluctuations and 
confirm the Yoshizawa hypothesis that when the turbu- 
lent field does not posses a unidirectional mean flow, the 
velocity turbulent transport term is not dominating the 
pressure transport (l2l . fl3l , H?| . In the present mixing 
both the advection and the production rate of the tur- 
bulent energy are zero and thus the turbulent transport 
(velocity and pressure) rate is of the same order of the 
dissipation rate. 



IV. INTERMITTENCY ASYMPTOTIC 
BEHAVIOR 

In this section we consider the asymptotic behavior 
with regards to the variation of the parameter that con- 
trols this kind of shearless mixing layer, that is the initial 
energy ratio £ = E\jEi between the high energy turbu- 
lent field 1 and the low energy turbulent field 2. As stated 
above, this ratio is unequivocally linked to the turbulent 
kinetic energy gradient. In this work, £ was varied be- 
tween 1.5 and 10 6 . The two external fields show, for 
moderate values of £, decay exponents which arc very 
close, so that the two homogeneous turbulences external 
to the mixing decay in a similar way and the value of 
E\jEi remains quite constant during the time interval 
considered H, @|. 

1 /2 

After few initial eddy turnover times t = i(Q)/E x (0) , 
where i is the initial integral scale (homogeneous through 
the whole domain) and -Bi(O) is the initial energy of the 
high energy side, a true mixing layer begins to emerge 
from the initial conditions and reaches a self-similar state. 
This means that all normalized moments distributions 
across the mixing collapse to a single curve when the 
position is normalized with the mixing layer thickness, 
which is defined as the distance between the points with 
normalized energy (E — E-2)I(E\ — E-i) equal to 1/4 and 
3/4, see sketch in fig.l. This definition has been used in 
many previous works on shearless mixing [2|, y, [5j . 

Results from numerical simulations show that the mix- 
ing layer is highly intermittent in the self-similar stage of 
decay, and its intcrmittency is dependant on £ . In or- 
der to analyze the flow intermittency, moments of the 
component u, that is the component in the direction 
of the flow of turbulent kinetic energy, were computed 
(the averages are computed by integrating over planes at 
x = const). A particular focus was placed on the skew- 
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FIG. 6: Temporal evolution of the maximum of the skewness 
— / — \ 3/2 

S — it 3 / ( it 3 J in the mixing for various energy ratios rang- 
ing from 12 to 10 4 . (a) Numerical experiments at Re\ = 45. 
Empty symbols refer to simulations in a (27r) 2 x 8tt domain, 
the others to simulations in a (2-k) 2 x Aty domain. The dashed 
line is the value of the reference skewness in a simulation of ho- 
mogeneous and isotropic turbulence carried out on the same 
computational domain, bars represent the maximum fluctu- 
ations of this skewness. (b) Laboratory data at Re\ = 44.5 
(perforated plate experiment) and Re\ — 78.1 (bar grid ex- 
periment) from wind tunnel experiments where a spatial de- 
cay is observed 0. The time in laboratory experiments has 
therefore been computed using Taylor's hypothesis, as t = 
d/U, where d is the distance from the grid and U is the mean 
velocity across the grids. 



ness S = u 3 / (u 2 ^ and kurtosis K = u 4 / (u 2 ^ . 

The velocity fluctuation u is responsible for the energy 
transport across the mixing. The skewness distribution is 
a principal indicator of intermittent behavior. It vanishes 
in homogeneous isotropic turbulent flows and thus it re- 
mains close to zero in the fields external to the mixing. 
The skewness takes a positive value within the mixing 
layer. Figure [6l^a) shows the time evolution of the max- 
imum of the skewness for four simulations with energy 
ratios between 12 and 10 4 . During the initial eddy turn- 
over times the skewness increases steadily, before bending 
at a time varying from 1.5 (£ = 12) to 4 {£ = 10 4 ). At 
this point the mixing layer enters a near self-similar stage 
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FIG. 7: Temporal evolution of the maximum of the kurtosis 
K = u 4 / (u 2 ^ in the mixing for various energy ratios rang- 
ing from 12 to 10 4 . (a) Numerical experiments at Re\ = 45, 
Empty symbols refer to simulations in a (2ir) 2 x 87r domain, 
the others to simulations in a (2n) 2 x 4-7T domain. The dashed 
line is the value of the reference kurtosis in a simulation of ho- 
mogeneous and isotropic turbulence carried out on the same 
computational domain, bars represent the maximum fluctu- 
ations of this kurtosis. (b) Laboratory data at Re\ = 44.5 
(perforated plate experiment) and Re\ = 78.1 (bar grid ex- 
periment) from wind tunnel experiments where a spatial de- 
cay is observed Q|. The time in laboratory experiments has 
therefore been computed using Taylor's hypothesis, as t — 
d/U, where d is the distance from the grid and U is the mean 
velocity across the grids. 



of evolution. Figure |^b) shows the time evolution of the 
maximum of the skewness in the V&W experiments, the 
3,3:1 perforated plate experiment, where £ = 6.27, and 
the 3 : 1 bar grid experiment, where £ = 6.19. Since in 
the laboratory all the statistics decay in space, we have 
estimated an equivalent temporal decay by using Tay- 
lor's hypothesis. The corresponding time in laboratory 
experiments has been computed as t = d/U, where d is 
the distance from the grid and U is the mean velocity 
across the grids @]. By comparing parts (a) and (b) of 
figure [5] one can see a good agreement. The distribution 
with the lowest value of £ in part (a), which is 12, start 
to bend at 1.5 - 1.7 eddy turn-over times and has val- 
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FIG. 8: (a) Maximum of the skewness and (b) estimate of 
the asymptotic kurtosis value as a function of the initial en- 
ergy ratio (the horizontal dashed line indicates the Gaussian 
reference value), (c) normalized position of the maximum of 
the skewness in the mixing layer as a function of the initial 
energy ratio. Note that one can expect a higher level of in- 
termittency in the data of 0] since these experiments had 
non-unity integral scale ratio (I1/I2 — 1.5). 



ues of S max approaching those of 0- Note that in the 
laboratory experiment the ratio of macroscales is about 
1.5 (this value is estimated by considering the finiteness 
of Re\ according to Sreenivasan (1998) [lg). This agrees 
with the finding 0, Q that if the gradient of kinetic en- 
ergy and macroscale are concurrent the mixing process 
is enhanced. In fact, one sees here that an higher energy 
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gradient, £ = 12, produces the same skewness than the 
gradient of scale associated with the lower energy gradi- 
ent, £ = 6.19, in the V&W experiment. In our numerical 
experiment, for the higher £ ratios, we note a sort of 
damped oscillation that appears beyond the first maxi- 
mum. This seems also to be shown by the 3 : 1 bar grid 
experiment, sec figure IHJb). 

The value of maximum skewness inside the mixing 
layer as a function of the energy ratio is depicted in figure 
ISJa). For values of E1/E2 lower than 10 2 it scales almost 
linearly with the logarithm of the energy ratio, which is 
in fair agreement with the scaling exponent of 0.29 found 
in @. 

Figure [7J shows the temporal evolution of the maxi- 
mum of the kurtosis inside the layer. Here again the 
comparison between our numerical data and the data of 
the V&W experiment is presented. The numerical and 
the laboratory results contrast well for comparable val- 
ues of Ei I Ei . A high peak is shown at the end of the 
formation time interval where the mixing process devel- 
ops. This peak is followed by a decrease, that could be 
interpreted as the fact that the more extreme intermit- 
tent turbulent events take place at the end of the forma- 
tion interval and before the self-similarity sets in. In the 
numerical experiments which last more time scale units 
than those in the laboratory, the decrease is followed by 
another damped increase-decrease cycle, as in the skew- 
ness case. The time asymptotic values were estimated by 
averaging over the last cycle. Note that data in figures 
[5] and [7] from laboratory experiments were obtained in 
the presence of concurrent gradients of integral scale and 
kinetic energy. Also in the kurtosis case, it can be ob- 
served that a higher energy gradient produces the same 
intermittency than a gradient of scale associated with a 
lower energy gradient [5], Q . 

The distribution of the peak of kurtosis inside the mix- 
ing is shown in figure [5Jb) . From this figure it can be 
noted that the kurtosis reaches very high values, much 
higher than the value of 3, that is the Gaussian reference 
value indicated in the figure by the dashed line. The kur- 
tosis asymptote is in fact close to 10.5, which indicates 
the presence in the mixing layer of extremely intense in- 
termittent events. 

A similar behavior of the skewness and kurtosis max- 
ima can be seen in the mixing penetration, defined as 
the instantaneous position along the x direction of the 
maximum of the skewness normalized with the instanta- 
neous mixing layer thickness A(t), see figure Htc). The 
penetration becomes constant in the self-similar evolu- 
tion. The penetration physically highlights the region of 
maximum intermittency, which is located in the low en- 
ergy side of the mixing layer. An increase of the energy 
ratio enhances the penetration of the high energy side 
into the low energy side. An asymptotic value of about 
1.2 A is obtained for £ —> 00, which gives an indication 
of the penetration of an isotropic turbulent field into a 
quiescent field. 

An alternative measure of the anisotropy is given by 



the velocity gradient statistics. We have computed the 
third and fourth order moments of both the longitudinal 
velocity derivative du/dx and transverse velocity deriva- 
tive du/dyi (no summation over i). These are so defined 

s du/dx = (du/dxf/{{du/dxYf/\ 
s du/dy , = {du/d Vt Y/({du/d yi yfi\ 1 = 1,2 

K du/dx = {du/dxY/{(du/dxf)\ 
K du/ayi = (du/dyi)*/((du/dyi) 2 ) 2 , i = 1, 2 

The averages are computed by integrating over planes 
at x = const. Figure [9] shows the time evolution of the 
peak of the longitudinal and transverse velocity deriva- 
tive skewness and kurtosis within the mixing. The figure 
includes, for comparison, the values measured in the two 
homogeneous and isotropic turbulent fields outside the 
mixing and the values deduced from figures 5 and 6 of 
the review by Sreenivasan and Antonia (1997) [l|| for 
Re x = 45. 

We observe that the temporal evolution of all these 
velocity derivative statistics during the mixing decay 
presents an initial transient which is very similar to that 
shown by the velocity statistics, the transient length is 
the same in the two cases and there is no lag. The max- 
imum values are always reached at t/r ~ 4 and increase 
with £. The longitudinal derivative moments are always 
larger than the transverse derivative moments, the dif- 
ference decreases with the increase of £. For instance, 
for £ = 10 4 , absolute values as high as 4 — 5 arc reached 
for the skewness, while values of 55 and 38 are measured 
for the longitudinal and transverse kurtosis, respectively. 
The anisotropy picture yielded by these velocity deriva- 
tive correponds to that of a a higher intermittency along 
the inhomogencous direction than across it. 



V. CONCLUSIONS 

We considered the simplest kind of turbulent shear- 
less mixing process which is due to the interaction of two 
isotropic turbulent fields with different kinetic energy but 
the same spectrum shape. This mixing is characterized 
by the absence of advection, production of turbulent ki- 
netic energy and an integral scale gradient. Such a situa- 
tion can be seen as the simplest form of turbulence inho- 
mogencity that can lead to a departure from Gaussianity. 
The study was carried out by means of Navier-Stokcs 
direct numerical simulations based on a fully dcaliascd 
Fourier-Galcrkin pscudospcctral method of integration. 
The data base was analyzed through single-point statis- 
tics involving the velocity and pressure fluctuations. 

We determined the temporal asymptotic behavior of 
the self-similar state. We also obtained the asymptotics 
for very high energy ratios between the isotropic turbu- 
lent fields which, through their interaction, initiate the 
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FIG. 9: (Color online) Temporal evolution of the peak of the 
velocity derivative statistics within the mixing, (a) Longitu- 
dinal velocity dervative skewness. (b) Longitudinal velocity 
derivative kurtosis. (c) Transverse velocity derivative skew- 
ness. (d) Transverse velocity derivative kurtosis. The in- 
stantaneous values in the high and low homogeneous regions 
external to the mixing are shown with the red (dashed) and 
green (dash-dotted) lines respectively. The violet (dotted) 
line is the reference value for Re\ — 45 deduced from figures 
5 and 6 in ref. 1181. 



mixing process. The infinite limit of the turbulent en- 
ergy ratio corresponds to the interaction of a region of 
isotropic turbulence with a relatively still fluid. In this 
limit the turbulent energy gradient reaches the maximum 
observable value associated to a given energy in the high 
energy side of the mixing. In this limit the mixing pene- 
tration is maximum and is as deep as 1.2 times the mixing 
thickness. 

We observed the intermittency and anisotropy of the 
mixings. Anisotropy was found to be mild for second 
order moments, on the contrary it was very intense in 
third and fourth order moments. The time asymptotic 
behavior of the anisotropy was almost independent of the 
turbulent energy ratio (i.e. turbulent energy gradient). 



The anisotropy observed through the third and fourth 
order moments of the velocity derivatives (longitudinal 
and transverse) is also very intense, but depends on the 
turbulent energy ratio. 

Despite having no gradient of integral scale, no mean 
shear and thus no advection and no production of turbu- 
lent kinetic energy, all mixings showed a departure from 
a Gaussian state for any turbulent energy ratio. This sig- 
nifies that the absence of these flow properties does not 
imply a condition of no intermittency. On the contrary 
the intermittency is highly dependant on the turbulent 
energy ratio between the two interacting fields. The in- 
termittency has a constant asymptote when this ratio 
approaches to infinity, which is consistent with the max- 
imum value of the turbulent energy gradient that can be 
asymptotically attained in this limit. It is deduced that 
the presence of a gradient of turbulent kinetic energy is 
a sufficient condition for the onset of intermittency. For 
any turbulent energy ratio we verified that the pressure 
transport is not negligible with regard to the velocity 
transport as in recirculating turbulent flows. 

In conclusion, by assuming that the interaction of two 
isotropic turbulent fields with different kinetic energy but 
the same integral scale is the non-homogeneous turbulent 
flow with the lowest level of dynamical complexity, we 
propose the hypothesis that the existence of a gradient 
of turbulent energy is the minimal requirement for Gaus- 
sian departure in turbulence, since there is experimental 
evidence that it is a sufficient condition to promote in- 
termittency 
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APPENDIX A: MEAN PRESSURE FIELD IN 
THE TURBULENT SHEARLESS MIXING FLOW 

The shearless turbulent mixing that we have studied 
is a flow where the average momentum is zero since the 
initial condition and the boundary conditions are such 
as to not generate a mean flow. It should be noted that 
the laboratory configuration, at least, those to date, is 
somehow different. In fact, when the two interacting 
homogeneous isotropic turbulent fields are generated by 
grids placed in a wind tunnel, a mean (homogeneous, i.e. 
shearless) flow is present in the normal direction to the 
mixing. However, it is true that an acceleration along the 
mixing direction could emerge if the initial gradients of 
mean pressure and turbulent kinetic energy do not com- 
pensate. As in the laboratory situation, this mean flow 
would remain homogeneous, thus also in this case the 
mixing would be shearless (i.e. devoid of the production 
of turbulent kinetic energy) . Let us first consider the av- 
eraged Navier-Stokes equations without the introduction 
of any model. The mean momentum equation is: 



d t U, + djUiUj = -(l/p)d l P - djWuJ + vV 2 Ui 



(Al) 



where the capital letters denote mean quantities, the 
small letters fluctuations, and the overline denotes the 
statistical average. For t = we have L 7 , = and the 
only non zero derivative is in the x direction, so that 
these equations reduce to 



d t U = -{l/p)d x P-d x \ 



(A2) 



where U is the mean velocity in the mixing direction. 
It can be seen that if the initial pressure gradient term 
balances the gradient of the part of the initial turbu- 
lent kinetic e_nergy associated to the fluctuations in the 
x direction (u 2 = 2/3K), the acceleration term dtU is 
zero. In such a situation, a mean field will be absent. 
On the contrary, for example in the hypothetical case of 
an initial kinetic energy gradient facing a zero pressure 
gradient, a mean homogeneous (without shear) flow will 
be generated. 

In the present numerical experiment the initial velocity 
field is first introduced. Then, as is standard practice, 
the code builds the pressure field by using the Poisson 
equation obtained from the divergence of the momentum 
balance. Periodicity conditions plus a condition fixing 
the average pressure po value in the entire domain are 
used. Since the field is incompressible, the divergence of 
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FIG. 10: (Color online) Part (a): profiles of the mean pres- 
sure and second order moment of the velocity component 
in the mixing direction, simulation with initial energy ratio 
E1/E2 = 60. Part (b): mean flow acceleration, gradient of 
the second order moment of the velocity component in the 
mixing direction and of the mean pressure. Simulation with 
initial energy ratio E1/E2 = 60. Air at standard conditions. 



Ui is zero and we obtain the following averaged equation: 



S7 2 P/p = -didjuju] - djUiUj 



(A3) 



At t = the fluctuating velocity field is statistically uni- 
form apart from in the x direction (note: it remains so 
during the mixing process). By also considering the sym- 
metries of the initial velocity field, and in particular the 
fact that, outside the mixing, the field is uniform, we 
obtain 



dL p /P = -d. 



2 TP, 



d x P/p = -d x u 2 



(A4) 



Consequently, by coming back to (|A2j) . one can see that 
no mean acceleration is generated at t = 0. Figure \W\ 
shows the terms in equation (A2) - the pressure and 
turbulent kinetic energy gradients and dtU - in two 
instants. We have considered the field configuration 
observed in the laboratory experiment by Veeravalli 
and Warhaft (1989, 3:1 perforated plate experiment, 
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air flow at standard conditions), which is actually the 
field configuration that we tried to reproduce in this 
numerical experiment. In particular, we have estimated 
the dimensional values of the pressure gradients and 
pressure difference between the high turbulent energy 
and low energy regions of the mixing. If (dP/dx) rnax = 
is the maximum value of the mean pressure gradient 
and AP = the pressure difference between the two 
homogeneous regions, we have at the initial instant of 
the simulations: 

Ei/Et = 6.6, (dP/dx) max = 1.39Pa/m, 



AP = 2.30 x 1CT2 Pa, 2A = 2cm 

Eh./ Ei = 40, (dP/dx) max = 1.60Pa/m, 
AP = 2.71 x 10~ 2 Pa, 2A = 2cm 

Ex/Ei = 60, (dP/dx) max = 1.62Pa/m, 
AP = 2.72 x 10~ 2 Pa, 2A = 2cm 

It can be observed that these pressure differences 
are very small. As a consequence, measurements in the 
laboratory should be very difficult. 



